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ABSTRACT 


The  n/m  Job-Shop  Problem  is  formulated  as  an  Integer  Program  and 
the  asymptotic  algorithm  of  integer  programming  is  applied.   A  simple 
way  of  getting  the  group  problem  is  obtained  and  some  results  of  the 
structure  of  the  group  problem  are  derived.   Two  problems  are  found  to 
arise  with  this  approach:   (1)  the  group  problem  is  larger  than  can  be 
handled  by  current  shortest  route  algorithms,  and  (2)  the  relaxed 
non-negativity  constraints  will  be  violated  by  at  least  one  optimal 
group  problem  solution.   Future  research  will  attempt  to  resolve  these 
problems. 


535.335 


Group  Theory  in  the  General  n/ni  Job-Shop  Problem 


Introduction 

Conway,  Maxwell  and  Miller  have  suininarized  past  research  on  the 

general  n/m  job-shop  problem  with  these  words  of  encouragement: 

"The  general  job-shop  problem  is  a  fascinating  challenge.   Although 
it  is  easy  to  state,  and  to  visualize  what  is  required,  it  is  extremely 
difficult  to  make  any  progress  whatever  toward  a  solution.  Many 
proficient  people  have  considered  the  problem,  and  all  have  come  away 
essentially  emtpy-handed.   Since  this  frustration  is  not  reported  in 
the  literature,  the  problem  continues  to  attract  investigators..."^ 

2 
In  this  paper,  the  group  theoretic  approach  to  integer  programming 

will  be  applied  to  this  problem.   We  will  formulate  this  problem  as  an 

Integer  Program  and  derive  a  group  problem  from  the  optimal  LP  solution 

ot  this  formulationo   This  problem  is  characterized  by  large  size  and 

considerable  structure.   For  example,  the  5/3  problem  given  in  Figure  1 

11 

yields  a  group  of  size  (40)  .   However,  the  structure  of  this  group  seems 

very  simple;  it  can  be  represented  by  11  equations  involving  24  variables 
whose  coefficients  are  all  0,  1  or  -1.   The  relevant  question  will  be  how 
to  use  the  structure  of  the  problem  to  allow  treatment  of  the  large  groups. 
A  second  problem  which  must  also  be  handled  is  that  solutions  to  the  group 


Conway,  R.W.,  W.L.  Maxwell,  and  L.W.  Miller,  Theory  of  Scheduling, 
Addison-Wesley,  1967,  p.  103. 

2 
Gomory,  R.E.,  "On  the  Relation  Between  Integer  and  Non-Integer  Solutions 

to  Linear  Programs,"  Proceeding  of  the  National  Academy  of  Sciences,  53, 

pp.  250-256,  (1965). 

Shapiro,  J.,  "Dynamic  Programming  Algorithms  for  the  Integer  Programming 
Problem — I:   The  Integer  Programming  Problem  Viewed  as  a  Knapsack  Type 
Problem,"  Operations  Research  1^,  pp.  103-121,  (1968),  and  "Group  Theoretic 
Algorithms  for  the  Integer  Programming  Problem  II:   Extension  to  a  General 
Algorithm,"  Operations  Research  16,  pp.  928-947,  (1968). 


In  section  1  of  this  paper  an  integer  programming  formulation  of 
the  Job  Shop  Scheduling  problem  is  presented.   We  then  relax  the  integrality 
requirements  in  this  formulation  to  obtain  an  LP  problem.   In  section  2 
we  show  that  the  solution  of  this  LP  problem  is  easily  obtained  by 
Critical  Path  Methods  and  discuss  which  variables  in  this  solution  will 
be  non-integer.   In  section  3  we  derive  a  Group  Problem  which  characterizes 
solutions  in  which  these  non- integer  variables  are  made  integer.   In 
section  4  we  discuss  the  structure  of  this  Group  Problem.   Section  5 
illustrates  our  results  with  a  numerical  example  and  section  6  contains 
some  concluding  remarks . 

1.   Integer  Programming  Formulation 

The  n/m  job-shop  problem  involves  n  jobs  which  must  each  be  processed 
on  m  different  machines.   The  sequence  in  which  a  given  job  is  processed 
on  each  machine  is  fixed  and  not  necessarily  the  same  for  2  different  jobs. 
For  example,  job  1  might  be  processed  first  on  machine  1,  then  2,  then  3, 
etc.,  while  job  2  is  processed  first  on  machine  3,  then  1,  then  2,  etc. 

The  processing  of  a  job  on  a  machine  is  called  a  task.   Tasks  have 
predetermined  process  times  (possibly  zero)  and  for  a  given  job,  all  pre- 
ceding tasks  must  be  complete  before  a  particular  task  can  be  started. 
The  problem  is  to  determine  a  sequence  for  processing  each  job  on  each 
machine  so  that  some  criterion  is  optimized.   Here  we'll  restrict  our 
attention  to  minimizing  make-span,  the  time  to  complete  all  jobs. 

A  given  problem  can  be  represented  by  an  m  x  n  matrix  of  task  process 
times  and  another  m  x  n  matrix  showing  the  processing  sequence  for  each 
job.   This  data  is  depicted  in  Figure  1  for  a  5/3  problem.   In  the  tables 


shown,  the  jobs  are  numbered  1  to  5  and  tasks  1  to  3  with  1  being  the 

first  task  in  time  performed  on  a  job.   Machines  are  labeled  A,  B,  and  C. 

The  second  table  relates  task  sequence  to  machines  by  giving  the  machine 

on  which  task  1  of  job  1  is  performed,  etc.   The  data  for  this  problem 

was  generated  randomly  and  will  be  used  throughout  the  paper  for  illustration. 


re  1 


Job/Task 
1 
2 
3 
4 
5 


1 
9 
3 
6 
5 
2 


2 
3 
3 
7 
5 
7 


3 
8 
0 
5 
5 
0 


Process  Times 


Job /Task 

1 

2 

3 

1 

B 

C 

A 

2 

A 

C 

B 

3 

B 

A 

C 

4 

B 

A 

C 

5 

B 

C 

A 

Seq 

uence 

There  have  been  several  linear  integer  programming  formulations  of  this 

3 
problem;  the  most  recent  and  compact  one  given  here  is  due  to  Wagner.   We 

let  t..  be  a  non-negative  variable  giving  the  start  time  of  the  j    task 

on  the  i   job  and  p   be  a  constant  process  time  for  the  same  task.  (p.. 

will  be  an  element  of  the  first  matrix  in  Figure  1.)   Even  though  t..  is  not 

required  to  be  integer  from  its  physical  interpretation  as  a  time,  we  will 

see  in  sections  2  and  3  that  because  of  the  problem  structure,  if  p..  is 

integer,  t..  will  always  be  integer  in  an  optimal  solution.   This  allows 

us  to  treat  t  .  as  though  it  were  required  to  be  integer  and  use  all-integer 


Wagner,  H.M.,  "An  Integer  Linear-Programming  Model  for  Machine  Scheduling," 
NRLQ,  6,  No.  2,  June  1959. 


programming  techniques.   Finally,  we  let  y   be  a  0-1  variable  equal  to 
1  if  job  i  preceeds  job  l   on  machine  k  or  0  if  £  preceeds  i  and  M  a 
sufficiently  large  positive  integer,  and  the  formulation  preceeds  as 
follows. 

First,  there  are  task  precedence  contraints  for  each  job. 

^1^   ^ij  ^  ^.j-1  -^  Pi.j-1  j+2....m.  i  =  l,...n 

Then,  constraints  which  sayj,  in  effect,  only  one  job  can  use  a  machine  at 


once. 


^ij  -  hp-^ip  '  A'^  ^   =  ^••••''" 


all  Ij  &  2,p  for  which  both 
tasks  are  processed  on 
machine  k, 

^2)  hp^  ^ij-Pij"  ^^-4>-^ 

In  these  constraints,  the  0-1  variable  and  the  large  positive  M  are  used 
in  the  familiar  way  of  representing  conjunctive  constraints  to  relax 

exactly  one  of  the  constraints  depending  on  which  job  is  first  on  machine  k 

k 
(and  hence  on  what  the  value  of  y   is). 

M  must  satisfy  the  following  relations,  where  T*  and  T^  are  the 

maximum  and  minimum  values  assumed  by  (t„   -  t..)  in  any  allowable  optimal 

•^   £p    ij 

solution. 

(2b)   M  >^  T*  +  P 

itp 

To  determine  T*  and  T^  we  must  have  an  upper  bound  on  the  objective  equation 
from  which  it  is  easy  to  deduce  an  allowable  range  of  optimal  values  for 
(t„   -  t..).   Note  that  although  we  have  assumed  one  constant  M  for 


notational  ease.  In  practice  we  could  actually  use  a  different  value  in 
each  equation  of  (2). 

Finally,  the  objective  can  be  represented  with  the  following  constraints; 


(3) 


min  z  =  t 


t^  >  t,   +  p,  i  =  1, .  <,  o  ,n 

f  —  im   *  im  * 


In  all,  there  will  be  (m  -  1) "n  equations  of  type  (1),  (n  -  l)n»m  of 
type  (2),  and  n  of  type  (3)o   There  will  be  -^ — 5 — ~*'"  ^"^  variables 
whose  values  represent  a  particular  sequence  of  jobs  on  each  machine  and 
n-m  non-negative  t.   variables  giving  the  start  time  of  each  task. 


2.   LP  Solutions  of  the  Problem 

In  this  section  we  will  show  what  form  the  optimal  LP  solution  to 
this  problem  will  take. 

Allowing  the  0-1  variables  y..  to  assume  non-integer  values  is 
equivalent  to  relaxing  (2).   In  this  case  there  is  a  simple  optimal  solution 
in  which 

^11  =  0 
(A)       t   =  t   j^  +  ^ii-l  J  ~  2,....in,  i  =  l,...,n 

t,  =  max  t. 
f        in 

All  slacks  in  (1)  will  be  zero,  and  one  slack  in  (3)  will  be  zero,  the 
rest  basic.   Figure  2  shows  a  network  interpretation  of  this  solution 
for  our  5/3  example,  while  Figure  3  depicts  the  machine  usage  implied  by 
this  solution.   Note  that  the  solution  is  clearly  not  feasible  for  the 


scheduling  problem  since  some  tasks  overlap  on  the  times  they  use  machines. 
There  are  a  number  of  alternate  optimum  values  for  the  0-1  y   variables 
in  (2)  and  for  the  slacks  associated  with  those  variables.   Consider  the 
general  form  of  one  of  these  equations  as: 


(5) 


tj^  =  Sj^  +  t2  +  P2  -  yM 


^2  =  ^2-^  ^-^Pl-  ^^-y>*^ 


where  we  have  added  slacks  s^  and  s-  and  simplified  notation  for  this 

k 
discussion  by  dropping  subscripts  on  y   and  changing  subscripts  on  t 

and  t„  .   t,  and  t-  are  task  start  times  whose  values  are  assumed  to  be 
£,p    1      2 

predetermined  by  (4),  p^  and  p-  are  constant  process  times,  and  M  is  a 
constant  as  previously  defined «   Since  the  slacks  s^  and  s.  and  the  0-1 
variable  y  appear  only  in  the  two  equations  shown  here,  we  can  set  values 
for  them  without  worrying  about  the  rest  o£  the  problem  equations «   Our 
objective  is  to  find  basic  values  to  these  variables  which  satisfy  (5) 
(and  hence  (2)).   Along  with  the  values  set  by  (4),  these  will  form  a 
basic,  optimal  LP  solution. 
One  possible  set  of  values  is: 

^  -^  ^2  -  ^1  -  Pi 

y _ 


(6)         s^  =  M  -  p^  -  p^ 


S2  =  0 

If  M  is  set  according  to  (2b)  then  y  will  be  positive,  one  of  the 
slacks  will  be  zero  and  the  other  non-negative  so  this  solution  is 


Figure  2 


20 


Graph  for  the  Optimal  LP  Solution  of  the  5/3  Example 
Nodes  correspond  to  tasks.   The  number  over  a  task  is  the  process 
time,  the  number  below  the  start  time.   The  number  beside  F  is  the 
value  of  the  objective  function » 
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basic  and  feasible. 

Note  that  this  set  of  values  implies  that  in  general  all  of  the 
0-1  variables  will  have  positive  fractional  values.   There  is  an  alternate 
set  of  optimal  values  in  which  many  of  the  0-1  variables  can  be  made 
nonbasic,  thus  reducing  the  size  of  the  group  problem  obtained.   These 

values  are: 

,  if  t  ^  t  +  p 


M 


1  -  ^2   *'2 
if  t^  1  t^  +  p^ 


,  otherwise 


^1  = 


WP2  '  ^f  h  -  ^2  -^  P2 

t^-t2-P2+M         .  if  t2  1  t^  +  p^ 

W~Pi~P2  »  otherwise 

'  «^2  -  h  -^  Pi 

,  otherwise 


It  will  always  be  possible  to  redefine  y  so  that  it  becomes  0  in  the 
case  when  it  would  be  1  and  we  will  assume  that  this  is  done.   It  can 
be  seen  that  these  values  are  also  basic  and  feasible. 

It  is  easy  to  see  what  we  are  doing  by  looking  at  Figure  3  for  our 
5/3  example.   There  will  be  one  0-1  variable  and  2  equations  like  (2) 
for  each  possible  pair  of  tasks  on  a  machine,  so  in  this  example  there 
are  22  0-1  variables.   The  values  in  (7)  correspond  to  making  y  0  or  1 
whenever  the  tasks  in  Figure  3  do  not  overlap  on  a  machine,  and  setting 
y  as  in  (6)  only  when  tasks  overlap.   Since  there  are  only  11  overlaps 


10 

in  our  example,  we  have  reduced  the  number  of  basic  y   's  by  half. 

As  we  shall  see,  this  corresponds  to  reducing  the  size  of  the  group 

22     11 
from  M   to  M   .   However,  we  do  loose  something  by  this  reduction, 

since  we  will  end  up  dropping  the  equations  of  type  (2)  corresponding 

It 
to  non-basic  y.-  's  as  part  of  dropping  non-negativity  on  basics  and 

hence  the  solution  to  the  group  problem  might  violate  some  of  these 

equations,  i.e.  produce  overlap  in  tasks  that  previously  didn't  overlap. 

In  this  event,  we  would  have  to  explore  other  solutions  of  the  group 

problem  to  find  one  which  is  feasible  in  all  equations.   Nevertheless, 

it  seems  reasonable  to  expect  that  this  will  occur  infrequently  and 

that  these  equations  are  likely  things  to  give  up  in  order  to  reduce 

the  size  of  the  group  problem.   Therefore,  we  will  assume  in  subsequent 

work  that  the  form  given  in  (7)  is  to  be  used  with  y..  *s  redefined  to 

make  those  with  the  value  1  nonbasic.   We  might  push  this  idea  still 

further  by  attempting  to  juggle  task  start  times  to  get  a  solution  with 

fewer  overlaps  and  no  larger  objective  value,  but  we  won't  bother  at  this 

time. 


3.   Derivation  of  the  Group  Problem 

In  this  section  we  will  follow  the  strategy  suggested  by  the  group 
theoretic  3pproach  to  integer  programming.  Beginning  with  the  optimal 
LP  solution  to  (1),  (2)  and  (3)  which  has  some  variables  at  non-integer 
values,  we  derive  a  group  problem  whose  solution  set  corresponds  to  the 
set  of  "corrections"  to  (1),  (2)  and  (3).  A  correction  is  a  set  of  non- 
negative  integer  values  for  non-basic  variables  which  force  all  basic 


4 
Gomory,  op.  cit.  Shapiro,  op.  cit. 
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variables  to  integer  values.   In  determining  a  correction  via  the 
group  problem  we  ignore  the  non-negativity  constraints  on  the  basic 
variables  so  the  solution  determined  may  still  be  infeasible.  The 
group  problem  is  determined  by  mapping  all  non-basic  columns  and  the 
right  hand  side  of  the  optimal  LP  tableau  into  their  fractional  parts. 

We  start  by  writing  (1),  (2),  and  (3)  without  the  objective  function 
in  partitioned  matrix  form  as  follows: 


(^^>   Vb^^21^R^Wf  ^^22^2R                     =^2 

W        (III)  83^X3  +  R3^T^  +  A3Y^      -IS3^      =  b3 

av)  Wb^^i'^r  +Vb        -^^b  =  ^ 

(^^    ^51^3-^  Sl^R  -^Vr           =N 


Where  (I)  represents  equations  (1),  (II)  equations  (3)  without  the 

objective,  (III)  equations  (2)  in  which  3.  y        is  basic,  (IV)  the  companions 

to  the  equations  (2)  represented  in  (III)  in  which  slacks  are  basic,  and 

V  equations  (2)  in  which  no  y   's  are  basic.   T  and  T  are  vectors  of  the 

jJc>  o       K 

basic  and  nonbasic  task  start  time  variables  T..,  T„  a  vector  including 

ij    F 

the  variable  t^  defined  in  (3)  and  the  n-1  basic  slacks  associated  with 

(3),  Y  and  Y  are  the  basic  and  non-basic  0-1  y.„  variables,  and  S._, 
D        K  Ix.  Id 

1  =  1,2,3,  and  S  _,  i  =  4,5,  represent  basic  and  non-basic  slacks. 

The  entries  of  matrices  B^  and  R . , ,  i  =  1,...,5,  B   ,  and  R   are 
all  0,  1  or  -1.   The  entries  of  A  ,  A  ,  and  A  are  all  M,  -M,  and  0. 
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The  entries  of  b  ,  1  =  1,...,5  will  be  process  times  or  process  times 

+  M  and  hence  are  all  integers „ 

By  the  way  we  define  y  in  (7)  and  by  (2)  we  have  A_  =  (M)*I  and 

A,  "  (-M)'I  where  •  indicates  scalar  multiplication.  Also,  note  that 

the  coefficients  of  t, ,  and  t.   in  once  equation  of  (2)  are  just  the 

negative  of  the  coefficients  in  the  companion  equation.   Thus,  if  we 

associate  corresponding  rows  of  8(111)  and  8(IV)  with  the  same  machine 

equations  we  will  have  B^,  =  -B, ,  and  R_,  =  -R, , „   Moreover,  assuming 

31     41      31     41  »        o 

m  >^  2  by  (3)  and  (4)  we  have  F   =  Oo 

Using  these  facts  the  optimal  LP  basis  for  (8)  can  be  written  as 
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and  it  is  easy  to  check  that  its  inverse  is 
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We  will  now  show  that  B   and  B^-  are  integer  matrices.   From  this 
it  can  be  seen  that  all  entries  in  (9) (I),  (II),  (IV),  (V)  are  integer 
so  that  when  we  take  their  fractional  parts  in  deriving  the  group 
problem  we  will  get  zeros o   Thus  they  can  be  dropped  from  consideration 
leaving  (9)  (III)  from  which  the  group  problem  will  be  obtained. 

First  note  from  (1)  that  if  we  write  the  equations  defined  by  (1) 

in  the  order  i  =  1,  j  =  2s.oos.mi,  i  =  2 ,  j  =  2,..05m,,..i  =  n^  j  =  2,...,m 

and  order  the  variables  t .  j,  as  ti,5.o.,t,  ,  t„,,...,t   then  B,,  will 

ij     11      In   21      nm       11 

have  the  form 


11 


0 


0 


0 


1 

0 

0 

0 

0 

1 

1 

0 

0 

0 

0 

-1 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

-1 

0 

0 


which  is  an  identity  matrix  with  a  few  -I's  just  below  the  diagonal. 
More  precisely, 

B^^  =  I  -  E 


where  E  =   e, .    and: 

ij  '  ' 


(10) 


^J  = 


0  or  1     j  =  i-1,  i  =  2 n 

0         otherwise 
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Lemitia  1:   Let  E  be  an  nxn  matrix  as  defined  by  (10).   Then  E  =  0, 

2       2 
Proof:   Let  E  »  ||e   ||.   Then,  using  (10) 

2     n 
^ij  '   j^fj^lK^jK  "  ^i,i-l^i-l,l-2       ^  "  ^'"■•»" 


so 

1  j  =  i-2,  i  =  3, . . . ,n 

2 

e   = 

0  otherwise 


p  0 

Now,  using  induction,  assume  E  =  ||e   ||  has  the  form 

.  1  j  =  i-Jl  ,  i  =  1-Hl,...,n 

(11)        eV  = 

0  otherwise 

and   snow  e,.      has   the  same   form. 


n 

£+1  £+1  I 

^j      \f/iK^Kj      ==   ^i,i=^lVl.i-£-l  ==  ^  =   2^^ "• 


1  j    =   i-£+l 

i  =   2+2.  ,.  , 

0  otherwise 


Hence  E     has    form   (11)    and  equals   0. 


Lemma   2:      B   j"  =    (I   +  E  +  E^  +    ...    e""''") 

Proof:  ^ii^iJ  =    (I   -  E)(I  +  E  +  E^  +   ...    e""""^) 

=  i-e  +  e-e2  +  e2+...  -e""^  +  e""^  -  e" 
=  I  -  e" 


=   I. 


Since  E  is  integer,  E-^  is  integer  and  B   is  integer, 


16 


Using  (3)  we  can  see  that  B„-  will  have  the  form 


22 


1 

0 

0 

0 

0 

1 

-1 

0 

0 

0 

1 

0 

-1 

0 

0 

1 

0 

0 

-1 

0 

-1 


where  the  first  column  corresponds  to  t,,  the  first  row  corresponds 
to  the  single  binding  equation,  and  subsequent  columns  correspond  to 
basic  slacks. 


Lemma  3:   B   =  I 


Proof:   Let  B^^  =  11^^^ I  I  and  8^2  =  l|b^j||.   Then 


b..  =  I  h.^\.   =  b-,b,  .  = 
ij    K=l     J       ^ 


1 
0 


J  =  1 

J  *   1 


b..  =  I   b.^b^, 
ij       iK  Kj 

i^l  ^-^ 


=  ^'il^lj  "^  ^i^lj 


1-1      J  =  1 

0  +  1       j  =  ii    1?^1 

0  +  0       j  i«  i.l 
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2  -1 

Again,  since  B   is  all  Integer,  B  .  will  be  also.   Thus  B   and 

B,_  are  integer  and  it  can  easily  be  seen  that  all  matrices  and  right 

hand  sides  in  equations  (9)  1,  11,  IV^  V  are  integer  and  they  can  be 

dropped  from  the  group  problem,  allowing  us  to  concentrate  on  (9)  III 

for  the  group  equations. 

Writing  III  as  a  group  problem,  after  multiplying  through  by  M  we 

get 

(12)         ^-^n\l\i-'Si^\'-hi\lhR  -  ^h^  =  -^i^^i  ^  ^3  ^^°^  ^> 

where  (mod  M-1)  means  each  row  is  taken  modulo  M,  and  T  ,  S^  ,  and  S 

are  integer. 

4.   Structure  of  the  Group  Problem 

In  this  section  we  will  look  in  greater  detail  at  the  structure  of 
the  group  problem  expressed  in  (12)  and  also  derive  the  linear  objective 
function  for  the  group  problem. 

A  single  equation  from  (8) (III)  may  be  written  as 

<^3>      ^^p  -  hi  -  4  -  ^^u  -  Pij "  M 

where  the  j    task  on  job  i  and  p   task  on  i    both  use  machine  k  and 
overlap  in  the  current  LP  solution.   S   is  the  slack  on  (2). 

In  the  optimal  LP  tableau  (9) (III) ,  (13)  will  be  basic  relative 

to  y.j  •   Hence  to  put  (13)  in  optimal  tableau  form  we  must  eliminate 

k 

all  other  basic  variables.   S   is  non-basic  by  the  way  we  constructed 

our  optimal  LP  solution  in  (7) .   If  j  =  1,  then  t..  =  t   =  0  by  (4) 
and  is  non-basic.   Otherwise,  in  general  t   /  0  and  is  basic. 
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In  this  case  we  write  a  subset  of  the  precedence  constraints  for 
job  i  given  in  (1)  as  follows: 


^±2  '   ''il  ■  ^11  Pil 


*^i3  '   ^±1  ■       '   ^12         "  Pi2 


^j  -\r^  '  'iJ-1  ^  ^i.1  "  ' 


where  we  have  introduced  S.  »  ^  -  l9  0..,j-l  as  slacks.   Adding  these 

IS" 


constraints  together  we  gets 

J-1      J-1 

(14)       t.^  -  t..  -  IS,  ^  I  p. 

ij    il     ,  ir     ,  ir 
•'  r=l      r-1 


Adding  (14)  and  (13)  yields 

J-^       k      k     J 

t   -  t^,  -   E  S^   "  S^„  -  My"'„  »   Z  p,   -  M 
pq    il     ,  ir    U          ^  U.  ,  ir 

"^  r=l  r=l 

Now,  if  q  =  1 ,  t   is  non-basic  and  we're  done.   If  not,  we  can 
repeat  the  same  process.   For  all  cases  we  can  write  our  result  in  the 
following  form: 

P-1            J-1  k      k  J  P-^   -M 

t„  ,  +   Z  S„   -  t,,  -   Z  S,   -  S^„  -  My^„  =  E  p.  -   E  p. 

11            .  Ir  il     ,lr  m          •'id     ,  ir  ^^Ir 

r=l           r=l  r=l  r=l 

j-1     p-1 
where  we  interpret   E   and   E   as  0  if  j  or  p  =  1 . 

r=l     r=l 

It 
All  variables  in  this  equation  are  non-basic  except  y.  •   If  we 

divide  by  M  we  get  the  optimal  LP  equation.   Taking  this  mod  1  gives  the 
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corresponding  equation  in  the  group  problem.   If  we  multiply  by  M  and 
take  the  result  mod  M  we  get  the  following  equivalent  group  equation: 

p-1  j-1  j       p-1 

r=l  r=l  r=l      r=l 

If  there  are  K  overlaps  of  jobs  on  machines  the  group  problem  will 
consist  of  K  of  these  equations  and  the  size  of  the  group  will  be  M  . 

We  will  obtain  the  objective  for  the  group  problem  is  a  similar 
manner.   Following  the  group  theoretic  approach  we  want  the  LP  objective 
expressed  in  terms  of  reduced  costs  and  non-basic  variables. 

One  of  the  constraints  in  (3)  will  be  binding,  say  for  i  =  b,  and 
correspond  to  the  longest  job.   In  this  equation  t^  is  basic. 

•"f  ■  ^fl  "  'bm  "  Pbn 

If  we  add  (14)  to  this  with  i  =  b  and  j  =  n  we  get 

j-1       n 

tr  -  S.,  -  t,  -  -   IS.   =  Z  p. 
f    fl    bl     ,  br     ,  br 
r=l      r=l 

or,  for  our  group  problem  objective 

n  j-1 

(16)  min  t.  =  l  p,      +  t,  T  +  I   S, 

f     1  br    bl     ,  br 
r=l  r=l 

where  we  dropped  S   since  it  doesn't  appear  in  (15).   Note  that  this 
objective  is  tantamount  to  minimizing  the  completion  time  of  the  job 
which  is  longest  in  the  LP  solution.   This  could  be  a  very  poor  objective 
to  use  if  another  job  became   longest  in  the  process  of  obtaining  an 
integerizing  correction.   This  occurrence  corresponds  to  one  of  the 
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basic  slacks  in  (3)  going  negative  and  illustrates  one  of  the  problems 
with  this  approach. 

(15)  and  (16)  give  us  considerable  insight  into  the  structure  of 
the  group  problem.  Note,  for  example,  that  all  coefficients  will  be 
0,  -1,  or  1.   Also,  notice  that  we  can  always  construct  an  optimal 
solution  to  the  group  problem  using  the  slacks  S   which  has  a  0  objective 
value  since  these  slacks  don't  appear  in  (16).   This  solution  will  always 
be  infeasible  in  terms  of  non-negativity  of  the  basic  variables  and 
Illustrates  another  problem  to  be  overcome  with  this  approach. 

We  now  have  a  constructive  procedure  for  easily  getting  the  group 
problem  by  inspection  from  a  CPM  solution  and  Gantt  chart  as  shown  in 
figures  2  and  3.   This  procedure  is: 

(1)  Identify  all  job  overlaps  on  machines. 

(2)  For  each  overlap,  write  an  equation  like  (15). 

(3)  Identify  the  longest  job  and  write  an  objective  equation  like 
(16). 

5.   Numerical  Example 

Using  the  procedure  just  outlined  we  can  derive  the  group  problem 
for  the  5/3  example  presented  in  figures  1,  2,  and  3. 

Step  (1).   Examining  figure  3  the  following  11  overlaps  can  be  identified: 

Machine  A  -  Jobs  3  and  4,  3  and  1 

Machine  B  -  5  and  4,  5  and  3,  5  and  1,  4  and  3,  4  and  1,  3  and  1 

Machine  C  -  5  and  2,  1  and  4,  4  and  3 
Step  (2).   The  group  equations  for  these  overlaps  are  written  below  in 

matrix  form. 
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Step  (3).   The  longest  job  in  the  LP  solution  is  job  1  and  from  this  we 
get  the  objective: 


min  t^  =  20  +  t^^  +  S^^  +  S^^ 


A  feasible  solution  with  t  =  35  is  known  for  this  problem  and 
using  this  value  it  can  be  determined  that  M  =  40  will  satisfy  (2b). 
This  implies  a  group  size  of  (40) 

6.   Conclusion 

When  the  group  theory  approach  is  applied  to  Job  Shop  Scheduling, 
the  group  problem  can  be  easily  obtained  by  inspection  from  the  CPM 
solution  of  the  associated  LP  problem.   The  group  problem  is  characterized 
by  considerable  structure  with  all  coefficient  0,  1,  or  -1. 

Two  problems  arise  when  we  apply  this  approach.   First,  the  size 
of  the  group  (M  )  is  generally  much  larger  than  can  be  handled  by 
present  shortest  route  algorithms.   Secondly,  since  we  can  construct 
at  least  one  optimal  solution  to  the  group  problem  which  violates  the 
relaxed  non-negativity  constraints  on  basic  variables,  it  appears  that 
the  relaxation  of  these  constraints  inherent  in  this  approach  will 
cause  difficulty. 

Future  research  in  this  area  will  address  itself  to  these  two  problems. 
The  special  structure  of  the  group  problem  will  be  studied  to  see  if  it 
can  be  exploited  in  developing  an  algorithm  for  solving  the  large  group 
problems.   In  connection  with  the  second  problem,  branch  and  bound 
approaches  will  be  investigated  for  searching  alternate  solutions  to 
the  group  problem  to  find  one  which  satisfies  basic  non-negativity. 


Shapiro,  op.  cit. 


23 

In  this  context  the  group  problem  could  be  used  for  getting  lower 
bounds  on  partial  solutions.   Relaxations  of  the  group  problem  which 
reduce  its  size  and  still  produce  a  strong  lower  bound  are  being 
investigated.   Finally,  the  use  of  Generalized  Lagrange  Multipliers 
in  branch  and  bound  procedures  will  be  investigated. 


Shapiro,  J.F.,  "Generalized  Lagrange  Multipliers  in  Integer  Programming," 
M.I.T.  Working  Paper  371-69,  March  7,  1969,  revised  January  5,  1970. 
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